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ABSTRACT 

We  consider  the  effect  Induced  on  the  data  processor  of  a  signal 
detection  system  when  the  underlying  noise  distribution  functions  are 
varied  about  their  nominal  values.  We  first  consider  the  detection  of  a 
time  varying  deterministic  signal  In  additive  noise,  and  then  extend  our 
results  to  a  more  general  situation  In  which  the  signal  possesses  a 
random  amplitude.  Our  results  characterize  a  class  of  contaminants  of  an 
arbitrary  nominal  distribution  over  which  the  data  processor  can  be 
designed  using  the  nominal  distribution,  and  It  Is  seen,  for  example, 
that  such  a  class  can  contain  distribution  functions  which  can  differ 
greatly  from  the  nominal  distribution. 

I .  INTRODUCTION 

Consider  the  discrete  time  detection  of  a  signal  In  corrupting  noise. 

It  is  well  known  that  a  typical  corresponding  detector  consists  of  a  data 
processor  followed  by  a  threshold  comparator.  In  fact.  It  Is  also  well 
known  that  under  a  variety  of  fidelity  criteria  the  data  processor  is 
characterized  In  terms  of  the  relevant  likelihood  ratio.  We  note  that  it 
is  often  convenient  In  practice  to  employ  an  appropriate  function  of  the 
likelihood  ratio  for  the  data  processor.  However,  in  order  to  present 
specific  results  about  the  general  situation  we  will  focus  attention  on 
the  likelihood  ratio  Itself. 

An  Important  factor  which  often  limits  the  employment  of  the 
likelihood  ratio  within  the  data  processor  Is  inexact  knowledge  of  the 
underlying  statistical  distributions.  One  way  of  resolving  this  difficulty 
Is  to  employ  a  robust  detector,  such  as  the  kind  obtained  from  Huber's 
results  (e.g.  see  [1]).  Numerous  authors  have  investigated  various  aspects 
of  robust  detection  [2-6].  While  robustness  approaches  can  achieve  the  -h.v 
desired  goal  of  desensitizing  the  data  processors  to  Inexact  statistical 
knowledge,  the  would-be  user  of  such  detectors  quickly  becomes  aware  of  £  | 
their  Inherent  limitations  (for  example,  the  extension  of  Huber's  results 
to  account  for  realistic  forms  of  dependency  and  effects  of  nonstatlonar-  V 
Ity,  the  question  of  the  loss  In  performance  occasioned  by  the  robustness,  • 
etc.).-  For  this  reason.  It  Is  beneficial  to  consider  situations  where  a 
simpler  approach  can  be  employed.  In  this  paper  we  consider  two  commonly  £ 
encountered  situations;  the  first  consists  of  the  detection  of  a  time  \  ° 

varying  deterministic  signal  In  additive  (not  necessarily  Gaussian)  noise,  £ 
whereas  the  second  generalizes  this  situation  to  Include  a  time  varying 
signal  with  random  amplitude.  Our  results  allow  determining  when  the 
likelihood  ratio  Is  Invariant  to  perturbations  In  the  underlying  distri¬ 
butions  from  their  nominal  values,  thus  admitting  the  employment  of  the 

Presented  at  the  Tuenty-Second  Annual  Allerton  Conference  on  Comrunioation, 
Control ,  and  Computing,  October  3-5,  1984 ;  to  be  published  in  the 
Proceedings  of  the  Conference.  85  10  i  i  i*fi 
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nominal  distributions  to  design  the  data  processor. 


II.  DEVELOPMENT 


The  situation  under  consideration  can  be  modeled  as  a  choice  between 
H0:  y0 
Ml:  yl 

where  p^  and  are  finite  measures  defined  on  the  Borel  sets  of  lRn.  Let 
X  be  a  o-flnlte  dominating  measure  for  Pg  and  p^,  e.g.  we  could  take  X  * 

\  (Pg  +  Pj).  Consider  first  the  detection  of  a  time  varying  deterministic 
signal;  that  Is,  for  a  fixed  element  6  of  IRn  we  set 


p| (B)  =  po(B-{0}) 

for  all  Borel  sets  B  of  lRn. 


Note  that  the  time  varying  signal  0  Is 


represented  as  a  point  In  IT.  We  define  the  real,  nonnegative,  Borel 

dyQ  dp, 

measurable  functions  fg  =  and  f^  *  ,  and  note  that  f^(»)  =  fg(--0) 

a.e.  with 
that  the 


c  lunwvivio  iq  -jj-  qiivi  i  ^  —  dX  *  a  u  llu mhm»  i^\w/  *0' 

respect  to  X.  We  also  let  S  *  {xe!Rn:  f0(x)f,(x)  >0},  and  note 
likelihood  ratio  A(*)  given  by  u 


A(x)  - 


f-|(x) 


exists  and  is  positive  for  all  xeS.  Note  that  S  contains  those  elements  of 

BRn  for  which  we  obtain  a  nondegenerate  test;  the  cases  where  A(x)  *0  or 
A(x)  *«  are  easily  checked. 

Let  IS  denote  the  class  of  all  real,  nonnegative,  Borel  measurable 

functions  defined  on  lRn;  we  will  call  an  element  of  «  a  density.  We  then 
note  that  a  wide  class  of  variations  in  the  nominal  density  fg(*)  is 

generated  by  considering  all  densities  f:  lRn  •♦•IR  such  that  f(*)  = 
(l-e)fg(*)  +eg(*)»  where  geW  and  e  Is  a  real  number  satisfying  0£e<l 

(cf.  the  e- contamination  class  of  [1]).  In  this  case  the  actual  density 
under  HQ  Is  not  necessarily  the  nominal  density  fQ(*)  but  rather  an 

arbitrary  element  of  ,  with  an  associated  ge«  and  ee[0,l).  In  the 
following  we  consider  only  those  g  belonging  to  r6.  We  thus  are  actually 
testing  between 


Hq;  ( 1  -e)  Pg+er0, 

Hj :  (1-e)  p^j  +  , 

where  .for  each  Borel  set  B  In  IRn, 

r0(B)  «  /  gdX  and  r] (B)  -  rQ(B-{0}). 

B 

Equivalently,  we  test  between  the  densities 


Hq:  (1-e)  fQ(*)  +eg(«) 

H-j :  (1-e)  fg(«-0)  +eg(*-0). 
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We  are  Interested  In  the  relevant  likelihood  ratio  A_  _(•)  given  by,  for 


v<*> s — 

Note  that  A„ 
9*£ 

that,  for  xeS, 


A  fgU-e)  +et9(x-8)-f0(x-6)] 


f0(x)  +e[g{x)-f0(x)] 

(x)  exists  for  all  xeS  via  our  restrictions  on  e.  Note  also 


8Aqe(x)  [fQ(x)  +e(g(x)-fQ(x))]  [g(x-e)-f0(x-e)] 

3e  [fQ(x) +e(g(x)-fQ(x))]2 

[fg(x-e)  +e(g(x-e)-f0(x-e))][g(x)-f0(x)] 
[f0(x)+e(s(x)-f0(x))]2 


We  therefore  have,  for  xeS, 


*  0  if  and  only  if 


g(x-e)  =  g(x) 
fo(x-0)  f0(x) 


This  is  true  if  and  only  if  there  exists  a  Borel  measurable  function  p:  IR 
-*•  IR  such  that  for  all  xeS  we  have  g(x)  =p(x)fQ(x)  and  p(x)  =  p(x-6).  Since 

the  above  equivalence  is  true  for  all  ee[0,l),  we  have  established  the 
following  result: 

Theorem  1:  Let  A  be  a  nonempty  open  subset  of  [0,1)  under  the  induced 
topology,  and  let  ge«.  We  then  have  A_  Jx)  *  A(x)  for  all  xeS  and  for 

9»e  n 

all  eeA  if  and  only  if  there  exists  a  Borel  measurable  function  p:  IR  -*IR 
such  that  for  all  xeS, 

g(x)  *  p(x)fQ(x)  and  p(x)  *  p(x-e). 

Note  that  the  above  result  characterizes  the  subclass  of  contaminants 

g( * )  for  which  the  likelihood  ratio  A„  (x)  is  invariant  for  each  xeS  as  e 

g,e 

ranges  over  nonempty  subsets  which  are  open  in  [0,1).  This  is  true 
regardless  of  the  nonempty  open  set  A  and  in  particular  we  can  simply  let  e 
be  any  number  in  [0,1).  The  data  processor  of  the  detector  could  then  be 
designed  based  on  the  nominal  density  with  no  degradation  in  performance 
for  such  contaminants  g. 

More  generally,  let  us  consider  the  extended  situation  where  is 
nominally  governed  by  the  (Borel  measurable)  density  f^ ( - )  given  by 

?i<x>  ■  /Ri.fo(x-t)d5(t)> 

where  for  a  fixed  subset  C  of  the  nonzero  Integers  the  finite  measure  y  is 

defined  on  the  Borel  sets  of  lRn  and  Is  purely  atomic  with  atoms 

{{c6}:c«C}.  Let  S  *  {xelRn:  fQ(x)f^(x)  >0}.  This  case  Includes  the 

previous  one  and  extends  it  to  the  situation  where  the  signal  Is  allowed 
to  possess  a  certain  form  of  random  amplitude.  As  before,  we  consider 
variations  In  the  nominal  density  fg(*)  which  lie  In  In  particular, 

let  getf  be  such  that  there  exists  a  Borel  measurable  function  p:  lRn-*IR 

with  g{x)  -  p(x)fQ(x)  and 


p(x)  *  p(x-e)  for  xeS.  We  then  note  that  the  likelihood  ratio  which 
corresponds  to  the  nominal  density  f0(»)  Is  given  on  5  by 

.  ^(x) 

A(x) '  f^y  • 

whereas  the  one  which  corresponds  to  the  variation  ?(•)  = (1 -e)f0( - )  +eg(*) 
Is  given,  for  all  xeS  and  for  all  ee[0,l),  by  u 

(x)  +  £  /  n  Cg(x-t)-f0(x-t)]dy(t) 


fQ(x)  +e[g(x)-f0(x)] 
fn(x-ce) +e[g(x-ce)-f*(x-ce)] 


ceC  fQ(x)  +e[g{x)-fQ(x)] 

for  some  countable  collection  of. real  constants  k£  which  depend  only  on  c, 
6,  and  p.  Note  that  It  then  follows  from  Theorem  1  that  for  each  xeS, 

A„  „<x)  *  A(x)  for  all  ce[0,l}.  We  have  therefore  established  the 
following  result: 

Theorem  2:  Suppose  ge«  and  that  there  exists  a  Borel  measurable  function 

p:  lRn  -*1R  with  g(x)  *  p(x)fn(x)  and  p(x)  =  p(x-e)  for  all  xeS.  It  then 
follows  that  u 

A„  (x)  *  A(x)  for  all  xeS  and  for  all  ee[0,l). 

The  above  result  Is  thus  an  extension  of  Theorem  1  to  a  more  general 
situation;  however,  we  retain  only  the  sufficiency  condition.  This,  of 
course.  Is  the  condition  we  often  might  want  to  employ  in  practice. 

Note  that  a  key  element  In  both  Theorem  1  and  Theorem  2  is  that  a 
contaminant  which  is  a  periodic  multiple  (with  period  equal  to  the  signal) 
of  the  nominal  density  leads  to  an  unperturbed  likelihood  ratio.  The 
results  of  Theorem  1  and  Theorem  2  provide  us  with  an  easily  used  tool  to 
construct  remarkedly  unintuitive  examples  of  noise  contaminants,  all  of 
which  employ  the  same  optimal  data  processor  for  a  wide  variety  of  fidelity 
criteria.  For  example,  let  n  *  1  and  0*0.1.  If 


f0(x)  *  i 


-{* 


if  xe[-l ,1] 

1 

otherwise 


l.e.  the  associated  random  variable  is  uniform  on  [-1,1],  and  If  the 
perturbation  of  fp  Is  given  by 


4.5x107  +  5x10'2  If  xe[-i  - 


for  1  = -10,-9,. . . ,10 


JO  If  xeC-co.-lJudH 

\  5xl0"2  otherwise 

then  the  resultant  likelihood  ratios  are  the  same  over  S*  [-0.9,1],  even 
though  the  two  densities  are  very  dissimilar.  Notice  in  this  example  that 
f(x)  Is  nowhere  equal  to  fg(x)  on  the  interval  [-1,1]. 

For  another  example,  let  n*l,  0*2,  and 


TM 


f0(»)  -  0.05  It.10>,0](*)  - 


C' 


-  5  - 

05  If  xe[-10,10] 
otherwise 


Let  D  be  a  Cantor  set  (see,  for  example,  [7])  In  [0,1]  with  Lebesgue 
measure  Ae(0,1).  If  the  perturbation  of  fQ  Is  given  by 


f(x) 


1 

25x 


1 


40- 20 A 


If  xeD  +  {21)  for  1=-5,-4,... ,4 
If  xe(-»,-10)u(10,»)  , 
otherwise 


then  the  resultant  likelihood  ratios  are  the  same  over  S  =  [-8,10],  even 
though  the  Cantor  set  results  In  a  perturbation  of  f  which  Is  obviously 
very  dissimilar  to  fQ.  Notice  that  the  maximum  of  the  perturbation  f  can 

be  made  to  exceed  any  preassigned  real  number. 


III.  CONCLUSION 


We  have  considered  the  effect  induced  on  the  data  processor  of  a 
signal  detection  system  when  the  underlying  distribution  functions  are 
varied  about  their  nominal  values.  Having  noted  the  crucial  role  played 
by  the  likelihood  ratio  in  the  detection  system,  we  have  therefore  focused 
attention  on  the  variation  in  the  likelihood  ratio  itself.  When  testing 
first  for  the  presence  of  a  time  varying  deterministic  signal,  and  then  in 
a  more  general  situation,  we  showed  that  the  likelihood  ratio  is  invariant 
over  a  class  of  perturbations  of  the  nominal  noise  distribution. 
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